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ABSTRACT 

Context: The joint X-ray and gamma-ray observations of GRB 090618 by a large number of satellites offer an unprecedented pos- 
sibility of testing crucial aspects of theoretical models. In particular, it allows us to test (a) in the process of gravitational collapse, 
the formation of an optically thick e'^e~ -baryon plasma self- accelerating to Lorentz factors in the range 200 < F < 3000; (b) its 
transparency condition with the emission of a component of 10^^"^^ baryons in the TeV region and (c) the collision of these baryons 
with the circumburst medium (CBM) clouds, characterized by dimensions of 10^^~^^ cm. In addition, these observations offer the 
possibility of testing a new understanding of the thermal and power-law components in the early phase of this GRB. 
Aims: We test the fireshell model of GRBs in one of the closest (z = 0.54) and most energetic (Eiso - 2.90 x 10^^ ergs) GRBs, namely 
GRB 090618. It was observed at ideal conditions by several satellites, namely Fermi, Swift, Konus-WIND, AGILE, RT-2 and Suzaku, 
as well as from on-ground optical observatories. 

Methods: We analyze the emission from GRB 090618 using several spectral models, with special attention to the thermal and power- 
law components. We determine the fundamental parameters of a canonical GRB within the context of the fireshell model, including 
the identification of the total energy of the e^e~ plasma, Efj , the Proper GRB (P-GRB), the baryon load, the density and structure 
of the CBM. 

Results: We find evidences of the existence of two different episodes in GRB 090618. The first episode lasts 50 s and is characterized 
by a spectrum consisting of thermal component, which evolves between kT - 54 keV and = 12 keV, and a power law with an 
average index y - 1.75 + 0.04. The second episode, which lasts for ~ 100s, behaves as a canonical long GRB with a Lorentz gamma 
factor at transparency of F = 495, a temperature at transparency of 29.22 keV and with characteristic size of the surrounding clouds 
of Rci ~ 10^^-16 cm and masses of ~ 10^2-24 g. 

Conclusions: We support the recently proposed two-component nature of GRB 090618, namely. Episode 1 and Episode 2, by using 
specific theoretical analysis. We further illustrate that the episode 1 cannot be considered to be either a GRB or a part of a GRB event, 
but it appears to be related to the progenitor of the collapsing bare core leading to the formation of the black hole which we call a 
"proto-black hole". Thus, for the first time, we are witnessing the process of formation of a black hole from the phases just preceding 
the gravitational collapse all the way up to the GRB emission. 

Key words. Gamma-ray burst: general — Gamma-ray burst: individual: GRB 090618 — Black hole physics 



1. Introduction 



After the disc overy of 
lites fclebesad el et all 



the GRBs by the Vela 



19731: [ Strong & Klebesadel 
5^^ 



satel- 

^- , 119741: 

IStrong et al.i. il974l: IStrong*. '1975*). the first systematic analy- 
sis on a large sample of GRBs was possible thanks to the ob- 
servations of the BATSE instrument on board the Compto n 
Gamma-Ray Observer (CGRO) satellite dMeesan et al.' 11992'). 
The 4BATSE catal og (iMeeganL 1 19971: |Paciesas et al, 1999; 
iKan eko et al.', "2006") consists of 2704 confirmed GRBs, and it 
is widely used by the science community as a reference for 
spectral and timing analysis on GRBs. One of the outcomes of 
this early analysis of GRBs led to the classification of GRBs 
as a function of their observed time duration. Tqq was defined 
as the time interval over which the 90% of the total BATSE 
background- subtracted counts are observed. The distribution of 
the Tqo duration was bi-modal: the GRBs with Tqo less than 2s 
were classified as "short" while the ones with T^q larger than 



2s were classified as "long" (Klebesadel", 19921; iDezalay et al.L 
Il992; Kouveliotou et al, 1993; Tavani, 1998). 

After the success of BATSE, a large number of space 
missions dedicated to the GRB observations were launched. 
Particularly significant was the discovery o f an additional pro - 
longed soft X-ray emission by Beppo-SAX (ICosta et al.Lll997l) . 
following the usual hard X-ray emission observed by BATSE. 
The Beppo-SAX observed emission was named as the "after- 
glow", while the BATSE one as the "prompt" radiation. The af- 
terglow allowed to pinpoint more accurately the GRB position in 
the sky and permitted the identification of their optical counter- 
part by space and ground based telescopes. The measurement of 
the cosmological redshift for GRBs became possible and their 
cosmological nature was firmly established (Ivan Paradijs et al.l 
.1997) . 

The Beppo-SAX and related results led to rule out liter- 
ally hu ndreds of theoretical models of GRBs, (see for a review 
.RuffiniL 1200 ll) . Among the handful of surviving models, there 



1 



L. Izzo etal.iGRB 090618 



was the one by iDamour & Ruffinil (Il975h based on the mass- 
energy formula of Black Holes. This model can naturally explain 
the energetics up to 10^"^"^^ erg, as requested by the cosmologi- 
cal nature of GRBs, through the creation of an -plasma by 
vacuum polarization pr ocesses in the Kerr-N ewman geometry 
(for a recent review see iRuffini et all l2010bl) . This model was 
proposed a few months after the presentation of the discovery 
of GRBs by Strong ( Strong, 1975) at the AAAS meeting in San 
Francisco. 

It soon became clear that, as suggested by Goodman and 
Paczynski (iGoodmanl 1 1 986l: IPaczy nski, 1986), the presence of a 
Lorentz gamma factor larger than 100 could overcome the prob- 
lem of opacity of the -plasma and justify th e y-ray emis - 
sion of GRBs at cosmological distances (see e.g. iPiranl l2005l) . 
That the dynamics of an -plasma with a baryon load with 
mass Mb would naturally lead to Lorentz gamma factor in th e 
range (10^ - 10^) was demonstrated bv IShemr& PiranI (Il990l) : 
iPiran et al.l (ll993h : lMeszaros etal.1 (Il993h . The general solution 
for a bar yon load B = Mbc^ l Efj between and 10"^ was ob- 
tained in "Ruffini et al." (^200u). The interaction between the ac- 
celerated baryons with the CBM, indicated by Meszaros & Rees 
(1 19931). w as advocated to explain the nature of the afterglow (see 
e.g. lPiranLflQQg and references therein). 

The unprecedented existence of such large Lorentz 
gamma factors led to formulate the Relativis tic Sp ace-Time 
Transformations paradigm for GRBs (Ruffini Tt al.L HoOlb). 
Such a paradigm made it a necessity to have a glo bal, instead of 
a piece wise, description of a GRB phenomenon dRuffini et al.L 
l2001bl) . This global description led to the conclusion that the 
emission by the accelerated baryons interacting with the CBM 
indeed occurs already in the prompt emission phase in a fully 
radiative regime. A new inter pretation of the burs t structure 
paradigm was then introduced (iRuffini et al.L 12001 ah : the exis- 
tence of a characteristic emission at the transparency of an e^e~- 
plasma, the Proper-GRB, followed by an extended- afterglow 
emission. The relative intensity of these two components is a 
function of the baryon load. It was proposed that the case of 
B < 10"^ corresponds to the short GRBs, while the case of 
5 > 3 X 10""^ corresponds to the long GRBs. 

This different parametrization of the prompt - afterglow ver- 
sus the one of the P-GRB - extended- afterglow could have orig- 
inated years of academic discussions. However a clear cut ob- 
servational evidence came from the Swift satellite, in favour of 
the second parametrization. The Norris-Bonnell sources, charac- 
terized by an initial short spike-like emission in the hard X-rays 
followed by a softer extended emission, had been indicated in the 
literature as short bursts. There is a clear evid ence that they be- 
long to a new cla ss of "disguised" short GRBs, (Bernardini et al.L 
120071: ICaito etal... 2009. 2010; de Barros et al.. 2011), where the 
initial spike is identified as the P-GRB while the prolonged soft 
emission occurring from the extended- afterglow emission in a 
CBM typically of the galatic halo. These sources have a baryon 
load 10""^ < B < 7x10""^: they are just long GRBs exploding in a 
particularly low density CBM of the order of 10"^ particles/cm^. 
This class of sources has given the first evidence of GRBs origi- 
nating from binary mergers, strongly suppo rted also from d irect 
optical observations (Bloom et al., 2006; Fo ng et af].l2010l) . 

It is interesting that, independent of the development of new 
missions, the BATSE data continue to attract full scientific inter- 
ests, even after the end of the mission in the 2000. Important 
inferences, based on the BATSE data, on the spectra of the 
early emiss i on of the GRB have been made by Ryde ( 2004) and 
iRyde et all (l2006h . They have convincingly demonstrated that 
the spectral feature composed by a blackbody and a power-law 



plays an important role in selected episodes in the early part 
of the GRB emission. They have also shown, in some cases, a 
power-law variation of the thermal component as a function of 
time, following a broken power-law behavior, see Fig.fTTl 

The arrival of the Fermi and other satellites allowed fur- 
ther progresses in the understanding of the GRB phenomenon 
in a much wider ene rgy range. Thanks to the Gamma-Ray 
Burst Monitor (GBM) ^M eegan et"aD. l2009 l) and the Large Area 
Telescope (LAT) ( Atwood ^et al.L l2009b . additional data are ob- 
tained in the 8 keV - 40 MeV and 100 MeV - 300 GeV energy 
range. It has allowed, among others, this first evidence of a GRB 
originating from the collapse of a core in the late evolu tion of a 
massive star, what we have cal led the Proto Black Hole (iRuffiniL 
120111: iPenacchioni etaD.l2012h . 

In the specific case of GRB 090618, it has been possible to 
obtain a complete temporal coverage of the emission in gamma 
and X-rays, due to the joint observations by Swift, Fermi, 
AGILE, RT-2/Coronas-PHOTON, Konus-WIND and Suzaku- 
WAM telescopes. A full coverage in the optical bands, up to 100 
days from the burst trigger, has been obtained. This has allowed 
the determination of the redshift, z = 0.54, of the source from 
spect roscopical identification of absorption lines dCenko et al.L 
12009 ) and a recent claim of a possible supernova emission -10 
days after the GRB trigger. This GRB lasts for ~ 150 s in hard 
X-rays, and it is characterized by four prominent pulses. In the 
soft X-rays there are observations up to 30 days from the burst 
trigger. 

We have pointed out in Ruffini etaP (1201 Oah that two diff'er- 
ent episodes are present in GRB 090618. We have also showed 
that while the second episode may fit a canonical GRB, the first 
episode is not expecte d to be eithe r a part of a GRB or an inde- 
pendent GRB (Ruffini eTaD.l201ll) . 

In the present paper we enter in the merit of the nature of 
these two episodes. In particular: 

- in Section 2, we describe the observations and data reduction 
and analysis. We obtain the Fermi GBM (8 keV - 1 MeV and 
260 keV - 40 MeV) flux fight curves, shown in Fig. [2l fol- 
lowing the standard data reduction procedure, and make a 
detailed spectral analysis of the main emission features, us- 
ing a Band and a power-law with an exponential high energy 
cut-ofl' spectral models. 

- in Section 3, after a discussion about the most quoted GRB 
model, the fireball, we recall the main features of the fireshell 
scenario, focusing on the reaching of transparency at the 
end of the initial optically thick phase, with the emission 
of the Proper-GRB (P-GRB). In Fig. [3] we give the theoret- 
ical evolution of the Lorentz F factor as a function of the 
radius, for selected values of the baryon load, corresponding 
to fixed values of the total energy E^^f . The identification 
of the P-GRB is crucial in determining the main fireshell 
parameters, which describe the canonical GRB emission. 
The P-GRB emission is indeed characterized by the tem- 
perature, the radius and the Lorentz F factor at the trans- 
parency, which are related with the E^^f^ energy and the 
baryon load, see Fig. |4l We then recall the theoretical treat- 
ment, the simulation of the light curve and spectrum of 
the extended-afterglow and, in particular, the determination 
of the equations of motion, the role of the E QuiTemporal 
Surfaces (EQTS) (Bianco_&Ruffiiiil 12004, 2005a), as wefl as 
the ansatz of the spect ral energy distr ibution in the fireshell 
comoving frame, (see iPatricelli et al.L [201 ll and references 
therein). 
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The temporal variability of a GRB light curve has been in- 
terpreted in some current models as due to internal shock 
(Rees & Meszaros, 1994). In the fireshell model instead such 
temporal variability is produced by the interaction of the 
ultra-relativistic baryons colliding with the inhomogeneities 
of the CircumBurst Medium (CBM). This allows to perform 
a tomography of the CBM medium around the location of 
the black hole formation, see Fig. \TU[ gaining important in- 
formation on its structure. These collisions are described by 
three parameters: the hcbm average density, the filling fac- 
tor "R, the dumpiness on scales of 10^^"^^ cm and average 
density contrast 10"^ < (Sn/n) < 10. We then refer also to 
the explanation of the observed hard-to-soft behavior due to 
the drop of the Lorentz T factor and the curvature eff'ect of 
the EQTS. We then recall the determination of the instanta- 
neous spectra and the simulations of the observed multi-band 
light curves in the chosen time interval, taking into account 
all the thousands of convolutions of comoving spectra over 
each EQTS leading to the observed spectrum. We also em- 
phasize how these simulations have to be performed together 
and optimized. 

in Section 4, we perform a spectral analysis of GRB 090618. 
We have divided the total GRB emission in 6 time inter- 
vals, see Table [T] each one identifying a si gnifica nt feature 
in the emission process, see also iRao et aP (l20lTI) . We have 
considered two diff'erent sp ectral models in the data fitting 
procedure: a Band model (IB and et al.L 11993) and on e by a 
blackbody plus a power-law component, following e.g. lRydd 
(l2QQ4l) . We find that the first 50s of emission are well-fitted 
by both models, equally the following 9s, from 50 to 59s. 
The remaining part, from 59 to 151 s, is fitted satisfactorily 
only by the Band model, see Table [T] 

in Section 5, we proceed to the analysis of GRB 090618 in 
the fireshell scenario. In Section 5.1, we attempt our first in- 
terpretation of GRB 090618 assuming it to be a single GRB. 
We recall that the blackbody is an expected feature in the the- 
ory of P-GRB. From the spectral analysis of the first 50 s, we 
find a spectral distribution consistent with a blackbody plus 
a power-law component. We have first attempted a fit of the 
source identifying these first 50s as th e P-GRB, see Fig. [6l 
We confirm the conclusion reached in iRuffini et"aD (I2010ali 
that this interpretation is not sustainable for three diff'erent 
reasons, based on: 1) the energetics of the source, 2) the time 
duration and 3) the theoretical expected temperature for the 
P-GRB. We then proceed, in the sub-section 5.2, to an in- 
terpretation of GRB 090618 as a m ulti component syst em, 
following the procedure outlined in iRuflSni et al.l (1201 ll) . in 
which we outline the possibility of the second episode be- 
tween 50 and 151s to be an independent GRB. 
We identify the P-GRB of this second episode, as the first 4s 
of emission. We find that the spectrum in this initial emission 
can be fitted by a blackbody plus a power-law component, 
see Fig. [HI Since this extra power-law component can be due 
to the early onset of the extended- afterglow, we take it into 
account to perform a fireshell simulation which is shown in 
Fig.O with an energy E^^f = 2.49 x 10^^ erg and a baryon 
load 5 = 1.98 ± 0.15 X 10"^ In Figs. llOll 11121 we report the 
results of our simulations, summarized in Tablel3] We notice, 
in particular, the prese nce of a strong ti me lag in this GRB. 
A detailed analysis, see lRao et al.l (l201ll) . about the time lags 
in the mean energy ranges of 35 keV, 68 keV and 125 keV, 
reports a quite large lag , ~ 7 s, in the first 50s of the emission 
which is unusual for GRBs, while in the following emission, 
from 51 to 151 s, the observed lags are quite normal, ~ 1 s. 



- in Section 6, we perform a spectral analysis of the first 50s, 
where we find a strong spectral variation with time, as re- 
ported in Table 13 and in Figs. I16I17I with a chacteristic 
power-law time var iation similar to the ones identified by 
iRvde&Pe'e'^ (12009 ) in a sample of 49 BATSE GRBs. 

- in Section 7, we estimate the variability of the radius emitter. 
Fig. HH and proceed to an estimate of the early expansion 
velocity. We interpret this data as originating in the expan- 
sion process occurring previous to the co llapse of the core 
of a m assive star to a black hole, see e.g. lArnett & MeakinI 
(12011 ): this early 50s of the emission are then defined as the 
proto-black hole phenomenon. 

- In Section 8, we proceed to the conclusions. 

2. Observations 

On 18th of June 2009, th e Burst Alert Telesc ope (BAT) on 
board the Swift satellite jGehrels et al. I l2009h triggered on 



GRB 090618 ( Schadv ltatl ^QQ9h. 



After 120 s 
and the 



the X-Ray 
Ultraviolet 



Telescope (XRT) (Burrows et al, 2005") 
Optical Telescope (UVOT) (Roming et all l2005h on board the 
same satellite, started the observations of the afterglow of 
GRB 090618. UVOT found a very bright opti cal coun t erpart 
with a white filter magnitude of 14.27 ± 0.01 (ISchadvl l2009l) 
not corrected for the extinction, at the coordinates RA(J2000) 
= 19:35:58.69 = 293.99456, DEC(J2000) = -^78:21:24.3 = 
78.35676. The BAT light curve shows a multi peak structure, 
whose total estimated duration is of ~ 320 s , with the Tqq dura- 
tion i n the (15-350) keV range was of 1 13 s (iBaumgartner et al.[ 
12009) . The first 50 s of the light curve presents a smooth decay 
trend, followed by a spiky emission, with three prominent peaks 
at 62, 80 and 112 seconds after the trigger time, respectively, and 
each have the typical appearance of the FRED pulse (see e.g. 
Fishman et al. (1994)), see Fig.l2j The time integrated spectrum, 
(to - 4.4, to + 213.6) s in the (15-150)keVrange, was found to be 
in agreement with a power-law spectral model with an exponen- 
tial cutoff', whose photon index w as y= 1.42 + 0.08 an d a cut-off* 
energy Epeak = 134 ± 19 keV (Sakamoto et al.V2009'). The XRT 
observatio ns started 125 s afte r th e BA T trigger time and lasted 
~ 25.6 ks (iBeardmore & Schadyil2009l) and reported an initially 
bright uncatalogued source, identified as the afterglow of GRB 
090618. Its early decay was very steep, ending at 310 s after the 
trigger time, when it starts a shallower phase, the plateau. Then 
the light curve breaks to a more steep last phase. 

GRB 090618 was observed also by the Ga mma-ray Burst 
Monitor (GBM) on board the Fermi satellite (Meegan et al.L 
I2009h . From a first analysis, the time-integrated spectrum, (^o, 
^0 + 140) s in th e (8-1000)keV range, was fitted by a Band 
(iBand et al.L 1 19931) spectral model, w ith a peak energ y Epeak = 
155.5 keV, a = -1.26 and/? = -2.50 (lMcBreenll2QQ9l) . but with 
strong spectral variations within the considered time interval. 

It is appropriate to compare and contrast the considerations 
of the time-integrated spectral analysis, often adopted in the cur- 
rent literature of GRBs, with the information from the time- 
resolved spectral analysis, as presented e.g. in this article. For 
a traditional astrophysical source, steady during the observa- 
tion time, the time-integrated and time-resolved spectral analysis 
usually coincide. In the case of GRBs, although the duration is 
only a few seconds, each instantaneous observation corresponds 
to a very diff'erent physical process and the two approaches have 
an extremely diff'erent physical and astrophysical content. 

The redshift of the source is z = 0.54 and it was determined 
thanks to the identification of the Mgll, Mg I and Fell absorption 
lines, using the KAST spectrograph mounted at the 3-m Shane 
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telescope at the Lick observatory (ICenko et al.ll2009l) . Given the 
redshift, and the distance of the source, we computed the emit- 
ted isotropic energy i n the 8 - 1000 keV energy range, using 
the Schaefer formula (ISchaefeii l2007l) : using the fluence in the 
(8-1000 keV) as observed by Fermi-GBM, S^^, = 2.7 x 10"^ 
(IMcBreenl l2009l) . and the ACDM cosmological standard model 
Hq = 70 km/s/Mpc, = 0.27, Qa = 0.73, we obtain for the 
isotropic energy emitted the value of E/^^ = 2.90 x 10^^ erg. 

This GRB was observed also by Konus - WIND 
(iGolenetskii et al.L l2009h . Suzaku-WAM (iKonoetaLl l2009l) 
and by the AGILE satelHte (iLongo et alll2009h . which detected 
emission in the (18-60) keV and in the MCAL instrument, 
operating at energies greater than 350 keV, but it did not observe 
high energy photons above 30 MeV. GRB 090618 was the first 
GRB observed by the Indian p ayloads RT-2 on board Russian 
Satell i te CORON AS -PHOTON (iKotov et al.Ll2008l:lNandi eTall 
120091: iRao et al.L 12011). Two detectors, namely, RT-2/S and 
RT-2/G consist of NaI(Tl)/CsI(Na) scintillators in phoswich 
assembly viewed by a photomultiplier tube (PMT). RT-2/S has 
a viewing angle of 4° x 4° and covers an energy range of 15 
keV to 1 MeV, whereas RT-2/G has an Al filter which sets 
the lower energy to ~ 20keV. The Mission was launched from 
Plesetsk Cosmodrom, Russia on January 30, 2009. During the 
event RT-2 payload was in the SHADOW mode (away from the 
Sun) during 08:16:10.207 UT and ended at 08:37:35.465 UT 
and the GRB 090618 was detected at 77° off'-axis angle. During 
this period, the spectra was accumulated in every 100s while 
the eight channel count rates for each detector are accumulated 
every second. The entire episode was observed for a duration 
of more than 200 seconds. A closer examination of the data 
in the accumulated Channels 1:15-102 keV, 2:95-250 keV and 
3:250-1000 keV indicates that the most significant counts is 
in Channel 2 with a clear evidence of the followings: (a) The 
emission in the first 50 s is prominent and broader in the lower 
channels, see Fig.[T] (b) After the first 50 s, there is an evidence 
of a precursor of about 6 seconds duration before the main pulse 
(c) a break up into two peaks of the main pulse at intermediate 
energies (35-200 keV) while at higher energies (250- 1000 keV) 
only t he first peak of the main pulse survives, see iRao et al.l 
(l201lh and also Fig. [2] in this manuscript. 

Thanks to the complete data coverage of the optical after- 
glow of GRB 090618, the possible presence of a supernova 
underlying the emission of the GRB 090618 optical afterglow 
(ICano et aH , 2011) was reported. The evidence of a supernova 
emission came from the presence of several bumps in the light 
curve and by the change in Rc - i color index over time: in the 
early phases, the blue color is dominant, typical of the GRB af- 
terglow, but then the color index increases, suggesting a pres- 
ence of a core-collapse SN. At late times, the contribution from 
the host galaxy was dominant. 



2.1. Data Analysis 

We consider the BAT and XRT data of the Swift satellite, to- 
gether with the Fermi-GBM and RT2 data of the Coronas- 
PHOTON satellite. The data reduction was done using the 
Heasoft v6.10 package^ for BAT and XRT, and the Fermi- 
Science tools for GBM. 

We obtained the BAT light curve and spectra using the stan- 
dard headas procedure. After the data download from the gsf c 
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Fig. 1. RT2 light curves of GRB 090618. 



websit^E we made a detector quality map and corrected the 
event data for the known errors of the detector and the hot pixels. 
We subtracted the background from the data, corrected for the 
improved position, using the tool batmaskwtevt and obtained 
the 1-s binned light curves and spectra in the main BAT energy 
band 15 - 150 keV and its subranges, using the tool batbinevt. 
After the systematic corrections to the spectrum, we created the 
response matrices and obtained the final spectra. 

For the XRT data, we obtained a total dataset using the 
standard pipeline, while for a time-resolved analysis we con- 
sider ed the on-l i ne recipe, w hich is well described in literature, 
see Evans et al.l (I2007ll2009 ). The GBM datfl in particular the 
fourth NaT detector in the (8 - 440 keV) and the bO BGO detector 
(260 keV - 40 MeV), were analyzed using the gtbindef tool to 
obtain a GTI file for the energy distribution and the gtbin for 
the light curves and final spectra. In order to obtain an energy 
flux lightcurve, we made a time resolved spectral analysis di- 
viding the count lightcurve in six time intervals, each of them 
corresponding to a particular pulse, as described in the work of 
iRao et al.l (l201ll) . All the time resolved spectra were fi tted us- 
ing the XSPEC data analysis software (lArnau version 
12.6.0q, included in the Heasoft data package, and considering 
for ea ch spectrum a classical Band spectral model (iBand et al.L 
fT993h and a power-law model with an exponential energy cut- 
off', folded through the detector response matrix. After the sub- 
traction of the background, we fit the spectrum by minimizing 
the between the spectral models described above and the ob- 
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Fig. 2. Fermi-GBM flux light curve of GRB 090618 referring to 
the Nal (8-440 keV, upper panel) and EGO (260 keV - 40 MeV, 
lower panel) detectors. 



served data, obtaining the best-fit spectral parameters and the 
respective model normalization. In Table [T] we give the results of 
our spectral analysis. The time reported in the first column corre- 
sponds to the time after the GBM trigger time ttrig = 267006508 
s, where the parameter w as not constrained, we used its aver- 
aged value, as delineated in lGuetta etaP (l201lh JS = -2.3 ± 0.10. 
We have considered the chi- square statistic for testing our data 
fitting procedure. The reduced chi- square = X^/^^ where N 
is the number of degrees of freedom (dof) which is A/^ = 82 for 
the Nal dataset and = 121 for the BGO one. 

For the last pulse of the second episode, the Band model is 
not very precise (^^ = 2.24), but a slightly better approximation 
is given by the power-law with an exponential cut-off', whose fit 
results are shown for the same intervals in the last two columns. 
From these values, we build the ffux light curves for both the 
detectors, which are shown in fig.[2l 



We turn now to the XRT which started to observe GRB 
090618 ~ 120 s after the BAT trigger. Its early data show a 
continued activity of the prompt emission, fading away ~ 200 
s after the BAT trigger time. Then the light curve is well approx- 
imated with a power-law decay. In view of the lack of soft X-ray 
data before the onset of the XRT, we cannot exclude a previous 
pulse in the X-ray light curve emission of GRB 090618. The fol- 
lowing shallow and late decay phases, well-known in literature 



(ISari et al.L 1 19991 : iNousek et all l2006h . will not be analyzed in 
this paper since we focus in the first 200 s of the GRB emission. 

3. A brief review of the fireshell and alternative 
models 

3.1. The GRB prompt emission in the fireball scenario 

A variety of models have been developed to theoretically explain 
the observational properties of GRBs. One of the most quoted is 
the fireball model (se e for a review Piran (2 005 )). The model 
was first proposed bylCavallo & ReesI (ll978h . lGoodmanl (Il986h 
and Paczynso (119861), who have shown that the sudden release 
of a large quantity of energy in a compact region can lead to 
an optically thick photon-lepton plasma and to the production 
of e'^e~ pairs. The total annihilation of the e'^e~ plasma was as- 
sumed, leading to a vast release of energy pushing on the GBM: 
the "fireball". 

An alternative approach, originating in the gravitational col- 
lapse to a black hole, is the fireshell model (see for a review 
iRuflSni et al. (2010^ and ([Ruffini, 201 1)). There the GRBs orig- 
inate from an optically thick electron-positron plasma in ther- 
mal equilibrium, having a total energy of E^^^. Such plasma is 
initially confined between the radius of a black hole rh and the 
dyadosphere radius 



rds = n 



2a 



^tot 



h/nieC 
rh 



1/4 



(1) 



where, a is the usual fine structure constant, fi and c the Planck 
constant and the speed of light, and the mass of the elec- 
tron. The lower limit of E^^^ coincides with Eiso- The condi- 
tion of thermal equilibrium assumed in this model as shown by 
Aksenov et al. ( 2007), diff'erentiates this approach from the alter- 
native ones (e.g. the one by iCavallo & Reeslll978i) . see Section 
3.2. 

In the fireball model, th e prompt emission, including the 
sharp luminosity variations (Ramirez-Ruiz & Fenimorei l2000h 
are du e to the prolonged and variable activi ty of the "inner en- 
gine" dRees & Meszaroslll994l:lPiranL[2QQ5b . The conversion of 
the fireball energy to radiation originates by shocks, either in- 
ternal (w hen faster moving matter takes over a slower moving 
shell, see iRees & M eszaros (1994)) or external (when the mov- 
ing matter is slo wed down by the external medium surround- 
ing the burst, see lRees & MeszarosI (Il992h ). Much attention has 
been given to the Synchrotron emission from relativistic elec- 
trons, possibly accompanied by SSC emission to explain the ob- 
served GRB spectrum. These processes were found to b e con- 
sistent with the obser vational data of many GRBs ( Tavanl [19961: 
iFrontera et al.l l2000h . However, several limitations have been 
evidenced in relation with the low energy spect ral slopes of 
time-integrated spectra (Crider et al., 1997; Preece etall 120021: 
iGhirlandaetaLl |2002, 2003; DaigneetaU 2009) and time- 
resolved spectra (iGhirlanda et al.. 20031) . Additional limitations 
on SSC ha ve also been poin ted out by iKumar & McMahonI 
(12008a') and 'Piran etaU (l2009h . 

The latest phases of the afterglow are described in 
the fireball model by assuming an equation of motion 
given by the Blandford-M c Kee s elf- similar power-law so- 
lution (iBlandford & McKeeL Il976h . The maximum Lorentz 
factor of the fireball is estimated from the temporal oc- 
currence of the peak of the optical emission, which is 
identified with the pea k of t he forward e xterna l shock 
emission (iMolinari et al.L 120071: iRvkoflf et al.L l2009h in the 
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Table 1. Time-resolved spectral analysis of GRB 090618. We have considered six time intervals, each one corresp onding to a 
partic ular emission feature in the light curve. We fit the GBM (8 keV - 10 MeV) observed emission with a Band model (iBand et al.L 
fl993h and a power-law function with an exponential cut-off'. In the columns 2-4 are listed the Band low energy index a, the high- 
energy 13 and the break energy E^^^^ , with the reduced chi-square value in the 6^^ column. In the last three columns are listed the 
power-law index y, the cut-off' energy E^^^ and the reduced chi-square value respectively, as obtained from the spectral fit with the 
cut-off' power-law spectral function. 

Time Interval a p ^^^^ (keV) Xband 7 (keV) 

~0^^50 -0.77!^-^^ -2.33!^-^^ VKU^ Til 0.91!^-^^ 180.9+^^-^ IXT 

50-57 -0.93+fj-48 -2.30 + 0.10 104.98+^f/ 1.22 Lll^fJ-g^ 168.3+^^1^ 1.22 

57 -68 -0-93 qI -2-43+^-^^ 264.0+JJ-^ 1.85 l-Org-g^ 340.5+J6-0 1.93 

68 -76 -1-05 -2.49+^-| 243.9+^-^ 1.88 311.0+38:6 1^90 

76- 103 -1-06 png -2.65+[;-J9 125.7+23-2V i^23 157.7+22:2 139 

103 - 150 -1-50 n 18 -2.30 + 0.10 101. ijn:^ 1.07 l-5Qo2o 102.8+^^-8 1.06 



thin shell approximation (ISari & PiranL [19991) . There have 
been developed partly alternative and/or complementary 
scenarios to the fireball model, e.g . the ones ba s ed on : 
quasi-thermal Comptonizatio n (Ghisellini & Celottj , 1 19991) , 
Com pton drag emission (IZdziarski et al.L Il99ll: IShemiL 
Il994i) . Synchrotron emission from a decaying magnetic 
field (iPe er & ZhangL [2006), jitter radiation (Medvedev, 
I2OOQI). Compton scatterin g of synchrotron se lf abs orbed photons 
(iPanaitescu & Me szarosl'2QQQ: Stern & Poutanen','2004), photo- 
spheric emission (Eichler & Levinson, 2000; Meszaros & Reed, 
'2000VMeszaros, 2002; Daigne & Mochkovitch, 2002; Gia nniol 
,2006; Rvde & Pe'er, 2009; Lazzati & Begelman, 2010). In par- 
ticular, Ryd e & Pe'erl (120091) pointed out that the photospheric 
emission overcomes some of the difficulties of pure non-thermal 
emission models. 



3.2. The fireshell scenario 

In the fireshell model, the rate equation for th e e'^e~ pairs and 
its dynamics have been given by Ruffini^tiD ([20001) (the pair- 
electromagnetic pulse or PEM pulse for short). This plasma en- 
gulfs the baryonic material left over in the process of gravita- 
tional collapse having mass Mg, still keeping thermal equilib- 
rium between electrons, positrons and baryons. The baryon load 
is measured by the dimen sionless parameter B = Mbc^ lEfj . It 
was shown dRuffiniL IT999 ) that no relativistic expansion of the 
plasma can be found for B > 10"^. The fireshell is still opti- 
cally thick and self-accelerates to ultrarelativistic velocities (the 
pair-ele ctromagnetic-baryonic pulse or PEMB pulse for short, 
iRuffiniL Il999 ). Then the fireshell becomes transparent and the 
Proper - GRB (P-GRB) is emitted ( Ruflani et all. 12001 ah . The fi- 
nal Lorentz gamma factor at transparency can vary in a vast 
range between 10^ and 10^ as a function of Efj and B, see 
Fig. [3 For the final determination it is necessary to integrate ex- 
plicitly the rate equation of the e^e~ annihilation process and 
evaluate, for a given black hole mass and a given e^e~ plasma 
radius , the reaching of the transparency condition iRuflSni et~aD 
([2QQQI) . see Fig.[l 

The fireshell scenario does not require any prolonged activ- 
ity of the inner engine. After transparency, the remaining ac- 
celerated baryonic matter still expands ballistically and starts to 
slow down by the collisions with the GBM, having average den- 
sity Hcbm- In the standard fireball scenario ( Meszaros, 2006), the 
spiky light curve is assumed to be caused by internal shocks. 
In the fireshell model the entire extended- afterglow emission 
is assumed to originate from an expanding thin shell enforc- 
ing energy and momentum conservation in the collision with 
the GBM. The condition of a fully radiative regime is assumed 



([Ruffini et al.L 12001a'). This, in turn, allows to estimate the char- 
acteristic inhomogeneities of the GBM, as well as its average 
value. 

It is appropriate to recall a further diff'erence between our 
treatment and the ones in the current literature. The com- 
plete analytic solution of the equations of motion of t he bary - 
onic shell has been developed (Bianco & Ruffini, 2QQ4j 200"5bb . 
while in the current literatu re usually the Blandford - McKee 
(iBlandford & McKed. L1976) self- similar solution has been un- 
critically adopted (e.g. Meszaro s^]etaD.[l 993: Sari. ll997L[l998 | : 
IWaxman, 1997; Rees & Meszaros, 1998; Granotetal., 19991; 
'Panait escu & Mes z aros , 1998; Gru zino v & Waxman, 1999; 
van Paradijs et al.l. |2Q0Q|: iMeszarosl |20Q2[) . The analogies and 
diff'erences betwe en the two approach es have been explicitly 
pointed out in Bia nco & Rufl5ni[ ([2QQ5a[) . 

From this general approach, a canonical GRB bolometric 
light curve composed of two diff'erent parts is deffned: the P- 
GRB and the extended- afterglow. The relative energetics of 
these two components, the observed temporal separation be- 
tween the corresponding peaks, is a function of the above three 
parameters Efj , B, and the average value of the ncbm\ the ffrst 
two parameters are inherent to the accelerator characterizing the 
GRB, i.e., the optically thick phase, while the third one is inher- 
ent to the GRB surrounding environment which gives rise to the 
extended-afterglow. If one goes to the observational properties 
of this model of a relativistic expanding shell, a crucial concept 
has been the introduction of the EQTS. In this topic, also, our 
model diff'ers from the ones in the literature for deriving the an- 
alytic expression of the EQTS from the analytic solutions of the 
equations of motion (Bianco & Ruffini, 2005a). 

In this paper, we assume Efj = Eiso- This assumption is 
based on the very accurate information we have on the luminos- 
ity and the spectral properties of the source. In other GRBs, we 
have assumed Efj > Eiso to take into account the observational 
limitations, due to detector thresholds, distance eff'ects and lack 
of data. 



3.3. The emission of the P-GRB 

The lower limit of Eff is given by the observed isotropic en- 
ergy emitted in the GRB, Eiso- The identiff cation of the energy 
of the afterglow and of the P-GRB determines the baryon load B 
and, from these, it is possible to determine, see Fig.[4l the value 
of the Lorentz F factor at transparency, the observed temperature 
as well as the temperature in the comoving frame and the labo- 
ratory radius at transparency. We can determine indeed from the 
spectral analysis of the P-GRB candidate, the temperature kTots 
and the energy emitted in the transparency Epqrb- The relation 
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Fireshell radius (cm) 



Fig. 3. The evolution of the Lorentz F factor until the trans- 



parency emission, for a GRB of a fixed 



= 1.22 X 10^5 



(upper panel), and Efj = 1 .44 x 10"^^, for diff'erent values of the 
baryon load B. This computation refers to a mass of the black 
hole of 10 Mq and a r = dr(ne± + n^^-)crT = 0.61, where ctt is 
the Thomson cross- section and the integration is over the thick- 
ness of the fireshell (lRuffiniLfT999h . 



between these parameters can not be expressed by an analytical 
formulation: they can be only obtained by a numerical integra- 
tion of the entire fireshell equations of motion. In practice we 
need to perform a trial and error procedure to find the set of val- 
ues which fit the observations. 

As we are going to see in the case of GRB 090618, the di- 
rect measure of the temperature of the thermal component at the 
transparency off'ers a very important new information in the de- 
termination of the GRB parameters. In the emission of the P- 
GRB two diff'erent phases are present: one corresponding to the 
emission of the photons when the transparency is reached, and 
the second is the early interaction of the ultra-relativistic protons 
and electrons with the CBM. A spectral energy distribution with 
a thermal component and a non-thermal one should be expected 
to occur. 



3.4. The extended-afterglow 

The majority of works in the current literature has addressed 
the analysis of the afterglow emission as due to various 
combinations of Synchrotron and Inverse Compton processes. 



(a) 

Do ppler blue-shif ted toward the observer 




(b) 



£^^^=1.44x10'^^ erg 
E*e°*= 1.68x10^^ erg 
Ee°*= 1 .77x1 0^^ erg 
E*e°*= 1.22x10^^ erg 




P-GRB 


Afterglow 


E^3°'^ 1.44x10^^ erg //Y 




. Ef= 1 .68x1 0^^ erg / / \ 




E*e°'= 1.77x10^^ erg / / / 









Fig. 4. The fireshell temperature in the comoving and observer 
frame and the laboratory radius at the transparency emission 
(panels (a) and (b)), the Lorentz F factor at the transparency 
(panel (c)) and the energy radiated in the P-GRB and in the af- 
terglow in units of Efj (panel (d)) as a function of the baryon 
load for 4 diff'erent values of Efj . 
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see e.g. iPiranl (l2Q05l) . It appears, how ever, that this descrip- 
tion is not fully satisfactory (see e.g. Ghirlan da et al.L 120031: 
iKumar & McMahon, 2008b; Piran et al, 2009). 

We have adopted in the fireshell model a pragmatic ap- 
proach by making the full use of the knowl edge of the equa- 
tions of motion, of the EQTS formulations (iBianco & Ruffinil . 
2005b) as well as of the correct relativistic transformations be- 
tween the comoving frame of the fireshell and the observer 
frame. These equations, that relate the four time variables, are 
necessary for the interpretation of the GRB data. They are: a) 
the comoving time, b) the laboratory time, c) the arrival time, 
and d) the arrival time at the detector corrected by the cos- 
mological eff'ects. This is the content of the Relative Space- 
Time Transformation s paradigm, es sential for the interpretation 
of GRBs data ( Ruffi ni et al.Ll2001bl) . Such a paradigm made it a 
necessity to have a gl obal, instead of a p iecewise, description of 
a GRB phenomenon (Ruffini et al., 2001b). This global descrip- 
tion led to a new inte rpretation of the burst structure paradigm 
(iRuflSni et al.L 12001 ah . As recalled in the Introduction, a new 
conclusion, arising from the burst structure paradigm, has been 
that the emission by the accelerated baryons interacting with the 
CBM is indeed occurring already in the prompt emission phase, 
just after the P-GRB emission. This is the extended- afterglow 
emission, which presents in its "light curve" a rising part, a peak, 
and a decaying tail. Following this paradigm, the prompt emis- 
sion phase is therefore composed by both the P-GRB emission 
and the peak of the extended-afterglow. 

To evaluate the extended-afterglow spectral properties, we 
have adopted an ansatz on the spectral properties of the emission 
in the collisions between the baryons and the CBM in the comov- 
ing frame. We have then evaluate all the observational properties 
in the observer frame by integrati ng on the EQT S. The initial 
ansatz of thermal spectrum (.RuflSni et al.L l2001ah . has been re- 
cently modified to 

^=f— V— f - (2) 

where a is a phenomenol ogical parameter define d in the comov- 
ing frame of the fireshell dPatricelli et al.Ll2QTl1) . determined by 
the optimization of the simulation of the observed data. It is 
well known that in the ultrarelativistic collision of protons and 
electrons with the CBM, collective processes of ultrarelativis- 
tic plasma physics are expected, no t yet fully explored and u n- 
derstood (e.g. Weibel instability, see lMedvedev & Loebl (Il999h ). 
Pr omising resu lts along t his line have been already obt ained 
bv lSpitkovskvl (l2008) and lMedvedev & Spitkovskvl (|200^ . and 
may lead to the understanding of the physycal origin of the a 
parameter in Eq. [2l 

In order to take into due account the filamentary, clumpy and 
porous structure of the CBM, we have introduced the additional 
parameter 91, which describes the fireshell surface filling factor. 
It is defined as the ratio between the effectiv e emitting area of the 
firesh ell Ae/f and its total visible area A^is (iRuffini et al.L l2002l 
I2005h . 

One of the main features of the GRB afterglow has been the 
observation of hard to soft spectral variation, which is generally 
absent in the first spike-like emission , which we have identifie d 
as the P-GRB.lBernardini et al.l(l2007h : rCaito et al.l(l2009Ll2010l) : 
Ide Barros et al.l (|201 ll) . An explanation of the hard-to-soft spec- 
tral variation has been advanced on the ground of two diff'erent 
contributions: the curvature effect and the intrinsic s pectral evo- 
lution. In particular, in the work oflPeng et al.l (1201 ll) the authors 
use the model developed in lOinl (l2002h for the spectral lag anal- 



ysis, taking into account an intrinsic Band model for the GRBs 
and a Gaussian profile for the GRB pulses, in order to take into 
account the angular effects, and they find that both causes pro- 
vide a very good explanation for the observed time lags. Within 
the fireshell model we can indeed explain a hard-to-soft spec- 
tral variation very naturally, in the extended- afterglow emission. 
Since the Lorentz F factor decreases with time, the observed ef- 
fective temperature of the fireshell will drop as the emission goes 
on, so the peak of the emission will occur at lower energies. This 
effect is amplified by the presence of the curvature effect, which 
has origin in the EQTS concept. Both these observed features are 
considered as the responsible for the time lag observed in GRBs. 

3.5. The simulation of a GRB light curve and spectra of the 
extended-afterglow 

The simulation of a GRB light curve and the respective spec- 
trum requires also the determination of the filling factor *R 
and of the CBM density ucbm- These extra parameters are 
extrinsic and they are just functions of the radial coordinate 
from the source. The parameter 91, in particular, determines 
the effective temperature in the comoving frame and the corre- 
sponding peak energy of the spectrum, while ricbm determines 
the temporal behavior of the light curve. It is found that the 
CBM is typically formed of "clumps" of width ~ 10^^"^^ cm 
and average density contrast 10"^ < (Sn/n) < 10 centered 
on the value of 4 particle s I cw? , see Fig. [TOl and clumps of 
masses Mdump ~ lo^^"^^ g Particularly important is the de- 
termination of the average value of ncbm- Values of the order 
of 0.1-10 particles/cm^ have been found for GRBs exploding 
inside star forming region galaxies, while values of the or- 
der of 10"^ particles/ cm^ have been found for GRBs explod- 
ing in galactic halos (iBernardini et al.L l2007l: ICaito et al.L 120091: 
de Barros et al., 20i3). The presence of such a clumpy medium, 
already predicted in pioneering works of Fermi in the theoret ical 
study of interstellar matter in our galaxy ( FermiL IT9491 1 1 9541) . is 
by now well-established both from the GRB observations and by 
additional astrophysical observations, see e.g. the circum-burst 
medium observed in novae (Sharaetal., 1997), or by theoret- 
ical considerations on supergiant, massive stars, clumpy wind 
(iDucci et al.L l2009l). Interesting are th e considerations by Arnett 
and Meakin (I Arnett & Meakinll201 ll) , who have shown how re- 
alistic 2D simulations of the late evolution of a core collapse 
show processes of violent emission of clouds: there the 2D simu- 
lations differ from the one in ID, which show a much more regu- 
lar and wind behavior around the collapsing core. Consequently, 
attention should be given also to instabilities prior to the lat- 
est phases of the evolution of the core, possibly giving origin 
to the cloud pattern observed in the CBM of GRB phenomenon 
(D. Arnett private communication). 

The determination of the ^ and ucbm parameters depends 
essentially on the reproduction of the shape of the extended- 
afterglow and of the respective spectral emission, in a fixed 
energy range. Clearly, the simulation of a source within the 
fireshell model is much more complex than simply fitting the 
N{E) spectrum with phenomenological analytic formulas for a 
finite temporal range of the data. It is a consistent picture, which 
has to fi nd the best value for the parameters of the source, the 
P-GRB dRuffini et al.L l2001al) . its spectrum, its temporal struc- 
ture, as well as its energetics. For each spike in the light curve 
are computed the parameters of the corresponding CBM clumps, 
taking into account all the thousands of convolutions of comov- 
ing spectra over each EQT S leading to the observed spectrum 
(iBianco & Ruffinil 120051)113) . It is clear that, since the EQTS 
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encompass emission processes occurring at different comoving 
times weighted by their Lorentz and Doppler factors, the "fit- 
ting" of a single spike is not only a function of the properties of 
the specific CBM clump but of the entire previous history of the 
source. Any mistake at any step of the simulation process aff'ects 
the entire evolution that follows and, conversely, at any step a 
fit must be made consistently with all the previous history: due 
to the non-linearity of the system and to the EQTS, any change 
in the simulation produces observable eff'ects up to a much later 
time. This brings to an extremely complex procedure by trial and 
error in the data simulation, in which the variation of the param- 
eters defining the source are further and further narrowed down, 
reaching very quickly the uniqueness. Of course, we cannot ex- 
pect the latest parts of the simulation to be very accurate, since 
some of the basic hypothesis on the equations of motion, and 
possible fragmentation of the shell, can aff'ect the procedure. 

In particular, the theoretical photon number spectrum to be 
compared with the observational data is obtained by an aver- 
aging procedure of instantaneous spectra. In turn, each instanta- 
neous spectrum is linked to the simulation of the observed multi- 
band light curves in the chosen time interval. Therefore, both the 
simulation of the spectrum and of the observed multiband light 
curves have to be performed together and simultaneously opti- 
mized. 

4. Spectral analysis of GRB 090618 

We proceed now to the detailed spectral analysis of GRB 
090618. We divide the emission in six time intervals, shown in 
Table[T] each one identifying a significant feature in the emission 
process. We then fit for each time interval the spectra by a Band 
model as well as by a blackbod y with an extra power-law com- 
ponent, following 'Rydd (l2004l) . In particular, we are interested 
in the estimation of the temperature kT and the observed energy 
flux ^obs of the blackbody component. The specific intensity of 
emission of a thermal spectrum at energy E in energy range dE 
into solid angle AD is 



Nonn. counts 



I{E)dE = 



-A^dE. 



(3) 



h^c^ exp(E/kT) - 1 

The source of radius R is seen within a solid angle AD = 
7iR^/D^, and its full luminosity is L = AnR^crT^. What we are fit- 
ting however is the background- subtracted photon spectra A{E), 
which is obtained by dividing the specific intensity I{E) by the 
energy E: 

k^L E^dE 



1(E) 

A(E)dE = -Y^dE = 



IcrikTfD^h^c^ cxp(E/kT) - 1 
l5(pobs E^dE 



(4) 



TT^ikTY exp(E/kT) - 1 ' 

where h, k and cr are respectively the Planck, the Boltzmann 
and the Stefan-Boltzmann constants, c is the speed of light and 
(pobs = LliAnD^) is the observed energy flux of the blackbody 
emitter. The great advantage of Eq. d?]) is that it is written in 
terms of the observables (pobs and T, so from a spectral fitting 
procedure we can obtain the values of these quantities for each 
time interval considered. In order to determine these parameters, 
we must perform an integration of the actual photon spectrum 
A{E) over the instrumental response R{i, E) of the detector which 
observe the source, where / denotes the diff'erent instrument en- 
ergy channels. The result is a predicted count spectrum 



Cn(0 : 
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A(E)R(i, E)dE, 
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Fig. 5. The two episode nature of GRB 090618. 



where Emin(i) and Emax(i) are the boundaries of the i-th energy 
channel of the instrument. Eq. (O must be compared with the 
observed data by a fit statistic. 

The main parameters obtained from the fitting procedure 
are shown in Table [2l We divide the entire GRB in two main 
episodes, as advanced in lRuflSni et aL lfeOll*): one lasting the first 
50 s and the other from 50 to 151 s after the GRB trigger time, 
see Fig. [5] It is easy to see that the first 50 s of emission, cor- 
responding to the first episode, are well fitted by a Band model 
as well as a black-body with an extra power-law model. Fig. (6] 
The same happens for the first 9 s of the second episode (from 
50 to 59 s after the trigger time). Fig. [71 For the subsequent three 
intervals corresponding to the main peaks in the light curve, the 
black-body plus a power-law model does not provide a satisfac- 
tory fit. Only the Band model fits the spectrum with good ac- 
curacy, with the exception of the first main spike (compare the 
values of in the table). We find also that the last peak can be 
fitted by a simple power-law model with a photon index y = 2.20 
±0.03, better than by a Band model. 

The result of this analysis points to a diff'erent emission 
mechanism in the first 50 s of GRB 090618 and in the follow- 
ing 9 s. A sequence of very large pulses follow, which spectral 
energy distribution is not attributable either to a blackbody or a 
blackbody and an extra power-law component. The evidence for 
the transition is well represented by the test of the data fitting, 
whose indicator is given by the changing of the;^^ (Ndof = 169) 
for the blackbody plus a power-law model for the diff'erent time 
intervals, see table [2l Although the Band spectral model is an 
empirical model without a clear physical origin, we do check its 
validity in all of the time-detailed spectra with the sole exception 
of the first main pulse of the second episode. The;^^ correspond- 
ing to the Band model for such a main pulse, although better than 
the one corresponding to the blackbody and power-law case, is 
unsatisfactory. We are now going to a direct application of the 
fireshell model in order to make more stringent the above con- 
clusions and reach a better understanding of the source. 



5. Analysis of GRB 090618 in the fireshell scenario: 
from a single GRB to a multi-component GRB 

5.1. Attempt for a single GRB scenario: the role of the first 
episode 

We first approach the analysis of GRB 090618 by assuming that 
we are in presence of a single GRB and attempt to identify its 
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Table 2. Time-resolved spectral analysis (8 keV - 10 MeV) of the second episode in GRB 090618. 





Time Interval (s) 


a 




EoikeV) 


^BAND 


kT(keV) 


r 


^BB+po 


A 


0-50 


-0.74 + 0.10 


-2.32 + 0.16 


118.99 + 21.71 


1.12 


32.07 + 1.85 


1.75 + 0.04 


1.21 


B 


50-59 


-1.07 + 0.06 


-3.18 + 0.97 


195.01 + 30.94 


1.23 


31.22 + 1.49 


1.78 + 0.03 


1.52 


C 


59-69 


-0.99 + 0.02 


-2.60 + 0.09 


321.74 + 14.60 


2.09 


47.29 + 0.68 


1.67 + 0.08 


7.05 


D 


69-78 


-1.04 + 0.03 


-2.42 + 0.06 


161.53 + 11.64 


1.55 


29.29 + 0.57 


1.78 + 0.01 


3.05 


E 


78 - 105 


-1.06 + 0.03 


-2.62 + 0.09 


124.51 +7.93 


1.20 


24.42 + 0.43 


1.86 + 0.01 


2.28 


F 


105-151 


-2.63 + -1 


-2.06 + 0.02 


unconstrained 


1.74 


16.24 + 0.84 


2.23 + 0.05 


1.15 



First episode - Band model 71™= -s^o! 4^.^ . First episode - BB + po model 




Fig. 6. Time-integrated spectra for the first episode (from to 50 s) of GRB 090618 fitted with the Band, = 1.12 (left) and 
blackbody + power-law (right) models, = 1.28. In the following we will consider the case of a blackbody + power-law model and 
infer some physical consequences. The corresponding considerations in the case of the Band model are currently being considered 
and will be published elsewhere. 



Second episc^de (50-59 s) - Band 




Second episode (50-59 s) - BB -i- po 




Fig. 7. Time-integrated spectra for the first 9 s of the second episode (from 50 to 59 s after the trigger time) of GRB 090618 fitted 
with the Band,;f^ = 1.23 (left) and blackbody + power-law (right) models, = 1.52. 



components in a canonical GRB scenario, based on the fireshell 
model. We first attempt the identification of the P-GRB emis- 
sion. We have already seen that the integrated first 50 s can be 
well-fitted with a black-body at a temperature kT = 32.07 ± 1.85 
keV and an extra power-law component with the photon index y 
= -1.75 ± 0.04, see panel A in Fig.[7]and TableO Being the pres- 
ence of a blackbody component the distinctive feature of the P- 
GRB, we have first attempted an interpretation of GRB 090618 
as a sin gle GRB with the first 50 s as the P-GRB Ruffini et al. 
(l2QlQah . We have first proceeded to evaluate if the energetics 
of the emission in the first 50 s can be interpreted as due to 
a P-GRB. The energy emitted by the sole blackbody is Ebb = 



8.35^Q3g X 10^^ ergs. Recalling that the isotropic energy of the 
entire GRB 090618 is Ei^o = (2.90 ± 0.02) x 10^^ ergs, we have 
that the blackbody component would be ~ 2.9 % of the total en- 
ergy emitted in the burst. This would imply, see lower panel in 
fig.m a baryon load B ~ 10"^ with a corresponding Lorentz F 
factor of ~ 800 and a temperature of ~ 52 keV. This value is in 
disagreement with the observed temperature kTobs = 32.07 keV. 

One may attempt to reconcile the value of the theoretically 
predicted GRB temperature with the observed one by increasing 
Ef^f^ . This would lead to an E^^^^ = 4 x 10^^ ergs and a corre- 
sponding baryon load of B ^ 10""^. This would imply three ma- 
jor discrepancies: a) there would be an unjustified complemen- 
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tary unobserved energy; b) in view of the value of the baryon 
load, and the corresponding Lorentz T factor, the duration of 
the extended-afterglow emission would be more t han an orde r 
of magnitude smaller than the observed 100 s (iBianco et all 
12008^) : c) the duration of this first 50 s is much longer than the 
one typically expect ed for all P-GRBs identified in other GRBs 
dRuflani etalj,l2007|), which is at maximum of the order of ~ 10 
s. We have therefore considered hopeless this approach and pro- 
ceeded to a diff'erent one looking for multiple components. 

5.2. The multi-component scenario: the second episode as 
an independent GRB 

5.2.1 . The identification of the P-GRB of the second episode 

We now proceed to the analysis of the data between 50 and 
150 s after the trigger time, as a canonical GRB in the fireshell 
scenario, namely the second episode, see Fig. \5\ (Ruffini et al., 
I2OIII) . We proceed to identify the P-GRB within the emission 
between 50 and 59 s, since we find a blackbody signature in 
this early second-episode emission. Considerations based on the 
time variability of the thermal component bring us to consider 
the first 4 s of such time interval as due to the P-GRB emis- 
sion. The corresponding spectrum (8-440 keV) is well fitted 
(X^ = 1.15) with a blackbody of a temperature = 29.22 + 2.21 
keV (norm = 3.51 ± 0.49), and an extra power-law compo- 
nent with photon index 7 = 1.85 ± 0.06, (norm = 46.25 ± 
10.21), see Fig. [8l The fit with the Band model is also accept- 
able (x^ = 1.25). The fit gives a low energy power-law index 
a = -1.22 ± 0.08, a high energy index j3 = -2.32 ± 0.21 and a 
break energy Eq = 193.2 ± 50.8, see Fig.[8l In view of the the- 
oretical understanding of the thermal component in the P-GRB, 
see Section 3.2, we shall focus in the following on the blackbody 
+ power-law spectral model. 

The isotropic energy of the second episode is Etso = (2.49 ± 
0.02) X 10^^ ergs. The simulation within the fireshell scenario is 
done assuming E^^f = Eiso- From the upper panel in Fig. |4] and 
the observed temperature, we can then derive the corresponding 
value of the baryon load. The observed temperature of the black- 
body component is = 29.22 ± 2.21, so that we can determine 
a value of the baryon load ofB= 1.98 ± 0.1 5x 10"^, and deduce 
the energy of the P-GRB as a fraction of the total . We so 
obtain a value of the P-GRB energy of 4.33+^-^^ x 10^^ erg. 

Now, from the second panel in Fig. |4] we can derive the ra- 
dius of the transparency condition, to occur at rtr = 1.46 x 10^^ 
cm. From the third panel we derive the bulk Lorentz factor of Tth 
= 495. We compare this value with the energy measured in the 
sole blackbody component of Ebb = 9.24^^^^ x 10^^ erg, and 
with the energy in the blackbody plus the power-law component 
of Ebb+po = 5.43^^-^^ X 10^^ erg, and verify that the theoretical 
value is in between these observed energies. We have found this 
result quite satisfactory: it represents the first attempt to relate 
the GRB properties to the details of the black hole responsible 
for the overall GRB energetics. The above theoretical estimates 
have been based on a non rotating black hole of 10 Mq, a to- 
tal energy of E^^f = 2.49 x 10^^ erg and a mean temperature 
of the initial plasma of e^e~ of 2.4 MeV, derived from the ex- 
pression of the dyadosphere radius, Eq.[TJ Any refinement of the 
direct comparison between theory and observations will have to 
address a variety of fundamental issues such as, for example: 1) 
the possible effect of rotation of the black hole, leading to a more 
complex dyadotorus structure; 2) a more detailed analysis of the 
transparency condition of the e'^e~ plasma, simply derived from 
the condition r = dr(ne± +n^^-)crT = 0.67 (iRuffiniLir999>) : 3) an 



photons keV 



s 




10 20 50 100 200 Energy(keV) 



Fig. 9. The fireshell simulation, green line, and the sole black- 
body emission, red line, of the time integrated (tO-h50, t0-h54 s) 
spectrum of the P-GRB emission. The sum of the two compo- 
nents, the blue line, is the total simulated emission in the first 4 
s of the second episode. 



analysis of the general relativistic, electrodynamical, strong in- 
teractions descriptions of t he gravitational collapse core leading 
to a black hole formation, fcherubini et al.L l2009l: iRuflSni et"all 
2003: .RuffinL. 19991 

5.2.2. The analysis of the extended-afterglow of the second 
episode 

The extended-afterglow starts at the above given radius of the 
transparency, with an initial value of the Lorentz F factor of Fq 
= 495. In order to simulate the extended- afterglow emission, we 
need to determine the radial distribution of the CBM around the 
burst site, which we assume for simplicity to be spherically sym- 
metric, we infer characteristic size of AR = 10^^"^^ cm. We have 
already recalled how the simulate of the spectra and of the ob- 
served multi band light curves have to be performed together and 
jointly optimized, leading to the determination of t he fundamen- 
tal parameters characterizing the CBM medium (iRuflSni et al.L 
2007). This radial distribution is shown in Fig.[TOl and is charac- 
terized by a mean value of < ^ > = 0.6 part/cm^ and an average 
density contrast with a < 6n/n > ^ 2, see Fig. [TOland Table 
m The data up to 8.5 x 10^^ cm are simulated with a value for 
the filling factor = 3 x 10"^, while the data from this value 
on with ^ = 9x10"^. From the radial distribution of the CBM 
density, and considering the 1/F effect on the fireshell visible 
area, we found that the CBM clumps originating the spikes in the 
extended- afterglow emission have masses of the order of IQ^^-^^ 
g. The value of the a parameter has been found to be -1.8 along 
the total duration of the GRB. 

In Figs. [TT]we show the simulated light curve (8-1000 keV) 
of GRB and the co rresponding spectrum, u si ng the spectral 
mode l described in (iBianco & RuffiniL |2004 . (iPatriceUi et al.L 

HoO). 

We focus our attention, in particular, on the structure of the 
first spikes. The comparison between the spectra of the first 
main spike (tO-h59, t0-h66 s) of the extended-afterglow of GRB 
090618, obtained with three different assumptions is shown in 
Fig. [T2I in the upper panel we show the fireshell simulation of 
the integrated spectrum (t0-h59, t0-h66 s) of the first main spike, 
in the middle panel we show the best fit with a blackbody and a 
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Fig. 8. On the left panel it is shown the time-integrated power spectra (8-440 keV) for the P-GRB emission episode (from 50 to 54 
s after the trigger time) of GRB 090618 fitted with the blackbody -i- power-law models, = 1.15, while on the right it is shown the 
fit with a Band model, = 1.25. 



power-law component model and in the lower panel the best fit 
using a simple power-law spectral model. 

We can see that the fit with the last two models is not satis- 
factory: the corresponding^^ is 7 for the blackbody -\- power-law 
and ~ 15 for the simple power-law. We cannot give the^^ of the 
fireshell simulation, since it is not represented by an explicit an- 
alytic fitting function, but it originates by a sequence of complex 
high non-linear procedure, summarized in Sec.[3l It is clear by a 
direct scrutiny that it correctly reproduces the low energy emis- 
sion, thanks in particular to the role of the a parameter, which 
was described previously. At higher energies, the theoretically 
predicted spectrum is aff'ected by the cut-off' induced by the ther- 
mal spectrum. The temporal variability of the first two spikes are 
well simulated. 

We are not able to accurately reproduce the last spikes of 
the light curve, since the equations of motion of the acceler- 
ated baryons become very complica ted after the first in terac- 
tions of the fireshell with the CBM (iRuffini et al.L l2007h . This 
happens for diff'erent reaso ns. First, a poss ible fragmentation of 
the fireshell can occur (Ruffini etali|2007|). Moreover, at larger 
distances from the progenitor the fireshell visible area becomes 
larger than the transverse dimension of a typical blob of matter, 
consequently a modification of the code for a three-dimensional 
description of the interstellar medium will be needed. This is un- 
like the early phases in the prompt emission, which is the main 
topic we address at the moment, where a spherically simmetric 
approximation applies. The fireshell visible area is smaller than 
the typical size o f the CB M clouds in the early phases of the 
prompt radiation, (llzzo et al., 2010). 

The second episode, lasting from 50 to 151 s, agrees with 
a canonical GRB in the fireshell scenario. Particularly relevant 
is the problematic of the P-GRB. It interfaces with the funda- 
mental physics issues, related to the physics of the gravitational 
collapse and the black hole formation. There is an interface be- 
tween the reaching of transparency of the P-GRB and the early 
part of the extended-afterglow. This connection has already been 
introduced in literature (Pe'er et al., 2010). We have studied this 
interface in the fireshell by analyzing the thermal emission at the 
transparency with the early interaction of the baryons with the 
CBM matter, see Fig.[9l 

We turn now to reach a better understanding of the meaning 
of the first episode, between and 50 s of the GRB emission. To 




1e+17 1.2e+17 
radius(cm) 

Fig. 10. Radial CBM density distribution in the case of GRB 
090618. The characteristic masses of each cloud are of the or- 
der of ~ io^^"^4 g and 10^^ cm in radii. 

Table 3. Final results of the simulation of GRB 090618 in the 
fireshell scenario 



Parameter 


Value 


^tot 


2.49 + 0.02 X 10^3 ergs 


B 


1.98 + 0.15 X 10-3 


To 


495 + 40 


kTth 


29.22 + 2.21 keV 


Ep-GRB,th 


4.33 + 0.28 X 10^1 ergs 


<n> 


0.6 partjcnv' 


< 6n/n > 


2 part/cnr' 



this end we examine the two episodes in respect to: 1) the Amati 
relation, 2) the hardness variation and 3) the observed time lag. 
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Table 4. Physical properties of the three clouds surrounding the 
burst site: the Distance from the burst site (2^^ column, the radius 
r of the cloud, 3^^ column, the particle density p, 4^^ column and 
the mass M in the last column 



Cloud 


Distance (cm) 


r (cm) 


p (#/cm-^) 


M(g) 


First 


4.0 X 10^^ 


1 X 10^^ 


1 


2.5 X 10^^ 


Second 


7.4 X 10^6 


5 X 10^^ 


1 


3.1 X 10^3 


Third 


1.1 X 10^^ 


2 X 10^^ 


4 


2.0 X 10^2 



fit of the light curve of GRB 090618 - 
li-GBM light curve counts data of GRB 090618 ^ 




theoretical simulation - 
Fermi GBM count spectrum 8-1000 keV ^ 




Energy (keV) 



Fig. 11. Simulated light curve and time integrated (tO+58, 
tO+150 s) spectrum (8-440 keV) of the extended-afterglow of 
GRB 090618. 



^ 0.1 
E 



count spectrum of GRB 090618 - 
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First spike (60-7S a) - power-law nmL»i eajiss^ 7EjeB . 
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6. The Amati relation, the HR and the time lag of the 
two episodes 

6.1. The first episode as an independent GRB? 

We first check if the t wo ep isodes fulfill separately the Amati 
relation, (lAmati et all l2002h . By using the Band spectrum we 
verify that the first episode presents an intrinsic peak energy 
value of Ep^ist = 223.01 ± 24.15 keV, while the second episode 
presents an Ep^2nd = 224.57 ±17.4 keV. The isotropic energies 
emitted in each single episode are Eiso,ist = 4.09 ± 0.07 x 10^^ 
ergs and Eiso,2nd = 2.49 ± 0.02 x 10^^ ergs, so we have that both 
episodes satisfy the Amati relation, see fig. \T3\ The fulfillment 
of the Amati relation of episode 2 was expected, being the sec- 
ond episode a canonical GRB. What we find surprising is the 
fulfillment of the Amati relation of the first episode. 



Fig. 12. Simulated time integrated (tO-h58, t0-h66 s) count spec- 
trum (8-440 keV) of the extended-afterglow of GRB 090618 
(upper panel), count spectrum (8 keV - 10 MeV) of the main 
pulse emission (tO-h58, t0-h66) and best fit with a blackbody -h 
power-law model (middle panel) and a simple power-law model 
(lower panel). 



We first examine the episode 1 as a single GRB. We notice 
a sharp rise in the luminosity in the first 6 s of emission. We 
therefore attempted a first interpretation by assuming the first 6 
s as the P-GRB component of this independent GRB, as opposed 
to the remaining 44 s as the extended- afterglow of this GRB. A 
value of the fit gives E'^^f = 3.87 x 10^^ ergs and 5 = 1.5 x lO'^. 
This would imply a very high value for the Lorentz fact or at the 
transparency of ~ 5000. In turn, this value would imply (iRuffiniL 
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1000 




Fig. 13. Position of the first and second component of GRB 
090618 in the Ep^i - Eiso plane respect the best fit of the 
Amati relation, as derived following the procedure described in 
('Capoz ziello & Izzol 12010 ). The red circle corresponds to the 
first emission while the green circle corresponds to the second 
one. 



Il999h a spectrum of the P-GRB peaking around ~ 300 keV, 
which is in contrast with the observed temperature of 58 keV. 
Alternatively, we have attempted a second simulation by assum- 
ing all the observed data be part of the extended-afterglow of a 
GRB, with a P-GRB below the detector threshold. Assuming in 
this case Etso = Efj , B = 10"^, and assuming for the P-GRB a 
duration smaller than 10 s, as confir med from the observations 
of all existing P-GRBs (iRuflani et al.l.l2007l) . we should obtain an 
energy of the P-GRB greater than 10"^ ergs/cm^s, which should 
have been easily detectable from Fermi and Swift. Also this sec- 
ond possibility is therefore not viable. We can then conclude 
generally that in no way we can interpret this episode either as 
a P-GRB of the second episode, as proved in paragraph 3.2 or, 
as proved here, as a separate GRB. We then conclude that the 
fulfillment of the Amati relation does not imply for the source to 
be necessarily a GRB. 

6.2. The HR variation and the time lag of the two episodes 

We finally address a further diff'erence between the two episodes, 
related to the Hardness-Ratio behavior (HR) and their observed 
time-lag. The first evidence of an evolution of the GRBs power- 
law slope indexe s with time was observed in the BATSE GRB 
photon spectra (Crider et al.L ll997h . In the context of the fireshell 
scenario, as recalled earlier, the spectral evolution comes out nat- 
urally from the evolution of the comoving temperature, the de- 
crease of the bulk Lorent z T factor and from the curvature eff'ect 
(iBianco & RuflSnil [20041) , with theoretically predicted values, in 
excellent agreement with observations in past GRBs. 

In order to build the HR ratio, we considered the data from 
three diff'erent instruments: Swift-BAT, Fermi-GBM and the 
CORONAS-PHOTON-RT-2. The plots obtained with these in- 
struments confirm the existence of a peculiar trend of the hard- 
ness behavior: in the first 50 s it is evident a monotonic hard- 
to-soft behavior, as due to the blackbody evolution of the first 
episode. For the second episode, the following 50 to 151 s of the 
emission, there is a soft-to-hard trend in the first 4 s of emission, 
and a hard-to-soft behavior modulated by the spiky emission in 
the following 100 s. For the HR ratio we considered the ratio 
of the count rate detected from a higher energy channel to that 




120 140 

Time (s) 




100 120 140 

Time (s) 



Fig. 14. Hardness-Ratio ratios for the Swift BAT data in two dif- 
ferent energy channels: HRl = cts(25-50 keV)/cts(15-25 keV), 
HR2 = cts(50-150keV)/cts(15-50keV). 



of a lower energy channel: HR = ctg(HE)/ctg(LE). In particular, 
we considered the count rate subtracted for the background, even 
when this choice provides bad HR data in time region dominated 
by the background, where the count rate can be zero or negative. 
For the Swift data, we consider the HR ratio for two diff'erent 
energy subranges: the HRl ratio shows the ratio of the (50-150 
keV) over the (15-50 keV) emission while the HR2 ratio shows 
the ratio of the (25-50 keV) over the (15-25 keV) emission, see 
Fig.m 

A similar trend is found for the Fermi-GBM Nal and RT-2 
instruments, see Fig {151 In particular, the HR from Fermi ob- 
servations was done considering the counts observed by the bO 
BGO detector in the range (260 keV - 40 MeV) and the ones ob- 
served by the n4 Nal detector in the range (8 - 260 keV). In Fig. 
[T5l it is shown the HR ratio for the Fermi observations, where 
we rebinned the counts in time intervals of 3 seconds. From this 
analysis we see that the HR ratio peaks at the beginning of each 
pulse, also for the second episode pulses, but each peak of the 
second episode pulses is softer than the previous one, suggesting 
that these pulses are consequential in the second episode and are 
in general agreement with the advance of a fireshell in the CBM. 
Since RT-2 data clearly show both the episodes up to 1 MeV it 
complements the results obtained by Swift (up to 200 keV) and 
FERMI (up to 440 keV) in the high and the most interesting en- 
ergy range. Hardness ratio plot of (250-1000 keV)/(8-250 keV) 
indicates that first phase of both episodes are the hardest. 
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Fig. 15. Hardness-Ratio ratio for the Fermi data. We considered 
the cts observed in the (260 keV - 40 MeV) energy range over 
the (8 - 260 keV) energy range. The time reported on the x-axis 
is in terms of the Mission Elapsed Time (MET). The presence of 
some negative data points is due to the presence of noise, in other 
terms the non-presence of GRB emission, in the background- 
subtracted BGO count light curve. 



Finally, the evident asymmetry of the first episode, supported 
by the observations of a large time lag in the high and low energy 
channels, see fig. [21 suggests a diff'erent process at work. There 
is a very signific ant softening of the first episode, as reported in 
lRaoetal.l(l2011h . where it is observed a large time lag between 
the 15-25 keV energy range and the 100-150 keV one: the high 
energy photons peak ~ 7 seconds before the photons detected in 
the 15-25 keV energy range. This large time lag is not observed 
in the second episode, where the lags are of the order of ~ 1 s. 

Motivated by these results, we proceed to a most accurate 
time-resolved spectral analysis of the first episode to identify its 
physical and astrophysical origin. 

7. A different emission process in the first episode 

7.1. The time resolved spectra and temperature variation 

One of the most significant outcome of the multi-ye ar work of 
Felix Ryde and his collaborators, (see e.g.lRyde et al.'(2010) and 
references therein), has been the identification and the detailed 
analysis of the thermal plus power-law features observed in a 
time limited intervals in selected BATSE GRBs. Similar fea- 
tures have been also observed recently in the data acquired by 
the Fermi satellite (Rvde et al.L 12010: .Guiriec et al.L ,2011) . We 
propose to divide these observations in two broad families. The 
first family presents a thermal plus power-law(s) feature, with a 
temperature changing in time following precise power-law be- 
havior. The second family is also characterized by a thermal 
plus power-law component, but with the blackbody emission 
generally varying without specific power-law behavior and on 
shorter time scales. It is our goal to study these features within 
the fireshell scenario, in order to possibly identify the underly- 
ing physical processes. We have already identified in Sec. 4 that 
the emission of the thermal plus power-law component char- 
acterizes the P-GRB emission. We have also emphasized that 
the P-GRB emission is the most relativistic regime occurring 
in GRBs, uniquely linked to the process of the black hole for- 
mation, see Sec. 5. This process appears to belong to the sec- 
ond family above considered. Our aim here is to see if the first 



Fig. 16. Evolution of the BB-Fpowerlaw spectral model in the 
yF(y) spectrum of the first emission of GRB 090618. It is ev- 
ident the cooling of the black-body and of the associated non- 
thermal component with the time. In this picture we have pre- 
ferred to plot just the fitting functions, in order to prevent some 
confusion. 



episode of GRB 090618 can lead to the identification of the 
above first family of events: the ones with temperature changing 
with time following a power-law behavior on time scales from 1 
to 50 s. We have already pointed out in the previous section that 
the hardness-ratio evolution and t he large time lag observed for 
the first episode (Rao et al.', '2011) points to a distinct origin for 
the first 50 s of emission, corresponding to the first episode. 

We have made a detailed time-resolved analysis of the first 
episode, considering diff'erent time bin durations in order to have 
a good statistic in the spectra and to take into account the sub- 
structures in the light curve. We have then used two diff'erent 
spectr al models to fit the observed data, a classical Band spec- 
trum (iBand et al.L Il993l) . and a blackbody with a power-law 
component. 

In order to have more accurate constraints on the spectral pa- 
rameters, we made a joint fit considering the observations from 
both the n4 Nal and the bO BGO detectors, covering in this way 
a wider energy range, from 8 keV to 40 MeV. To avoid some bias 
due to low photon statistic, we considered an energy upper limit 
of the value of 10 MeV. We report in the last three columns of the 
Table [5] also the spectral analysis performed in the energy range 
of the BATSE LAD instrument (20-1900 keV), as analyzed in 
iRyde & Pe'ed (l2009l) . just as a comparison tool with the results 
described in that paper. Our analysis has been summarized in 
Figs. O [T7] and in Table [5l where we report the residual ratio 
diagram as well as the reduced-;^^ values for the spectral models 
considered. 

We conclude that both the Band and the proposed blackbody 
-h power-law spectral models fit very well the observed data. 
Particularly interesting is the clear evolution in the time-resolved 
spectra, corresponding to the blackbody and power-law compo- 
nent, see Fig.[T6j In particular the kT parameter of the blackbody 
presents a strong decay, with a temporal behavior well described 
by a double broken power-law function, see upper panel in Fig. 
[TTl From a fitting procedure we obtain the best fit (R^- statistic 
= 0.992) for the two decay indexes for the temperat ure variation 
are aki = -0.33 ± 0.07 and bkr = -0.57 ± 0.11. In Ryde & Peld 
( 2009) an average value for these parameters on a set of 49 GRBs 
is given: (akr) = -0.07 ± 0.19 and (bkr) = -0.68 + 0.24. We note 
however that in the sample considered in iRyde &Pe'ed (120091) 
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Table 5. Time-resolved spectral analysis of the first episode in GRB 090618. We have considered seven time intervals, as described 
in the text, and we used two spectral models, whose best-fit parameters are shown here. The last three columns, marked with a LAD 
subscript, report the same anal ysis but in the energy range 20 - 1900 keV, which is the same energy range of the BATSE-LAD 
detector as used in the work of lRyde & Pe'ej (120091) . 



Time 




a 




^0 (keV) 


^BAND 


kT (keV) 


r 


^BB+po 


kTiAD (keV) 


yiAD 


^BB+po,LAD 


A:0- 


5 


-0.45 ±0.11 


-2.89 ± 0.78 


208.9 ±36.13 


0.93 


59.86 ± 2.72 


1.62 ±0.07 


1.07 


52.52 ± 23.63 


1.42 ±0.06 


0.93 


B:5- 


10 


-0.16 + 0.17 


-2.34 ±0.18 


89.84 ± 17.69 


1.14 


37.57 ± 1.76 


1.56 ±0.05 


1.36 


37.39 ± 2.46 


1.55 ±0.06 


1.27 


C:10 


- 17 


-0.74 ± 0.08 


-3.36 ± 1.34 


149.7 ±21.1 


0.98 


34.90 ± 1.63 


1.72 ±0.05 


1.20 


36.89 ± 2.40 


1.75 ±0.06 


1.10 


D:17 


-23 


-0.51 ±0.17 


-2.56 ± 0.26 


75.57 ± 16.35 


1.11 


25.47 ± 1.38 


1.75 ±0.06 


1.19 


25.70 ± 1.76 


1.75 ±0.08 


1.19 


E:23 


-31 


-0.93 ±0.13 


unconstr. 


104.7 ±21.29 


1.08 


23.75 ± 1.68 


1.93 ±0.10 


1.13 


24.45 ± 2.24 


1.95 ±0.12 


1.31 


F:31 


-39 


-1.27 ±0.28 


-3.20 ± 1.00 


113.28 ±64.7 


1.17 


18.44 ± 1.46 


2.77 ± 0.83 


1.10 


18.69 ± 1.89 


4.69 ± 4.2 


1.08 


G:39 


-49 


-3.62 ± 1.00 


-2.19 ±0.17 


57.48 ± 50.0 


1.15 


14.03 ± 2.35 


3.20 ± 1.38 


1.10 


14.71 ± 3.52 


3.06 ± 3.50 


1.09 



only few bursts shows a break time around 10 s, as it is in our 
case, see Fig. [T71 There are two of these bursts, whose analy- 
sis presents many similarities with our presented source GRB 
090618: GRB 930214 and GRB 990102. These bursts are char- 
acterized by a simple FRED pulse, whose total duration is ~ 40 
s, quite close to the one corresponding to the first episode of 
GRB 090618. The break time tt in these two burst are respec- 
tively at 12.9 and 8.1 s, while the decay indexes are akr = -0.25 
± 0.02 and bki = -0.78 ± 0.04 for GRB 930214 and akr = -0.36 
+ 0.03 and bur = -0.64 ± 0.04 for GRB 990102, see Table 1 in 
iRvde&Pe'erl (l2009l) . in very good agreement with the values 
observed for the first episode of GRB 090618. We conclude that 
the values we observe in GRB 090618 are very close to the val- 
ues of these two bursts. We shall return to comp are and contrast 
our results with the other sour ces considered in (iRyde & Pe'ed. 
2009), as well as GRB 970828 (iPe'er etal.Lr2007h in a forthcom- 
ing publication. 

The results presented in Figs. 1161171 as well as in Tabled 
point to a rapid cooling of the thermal emission with time of 
the first episode. The evolution of the corresponding power-law 
spectral component, also, appears to be strictly related to the 
change of the temperature kT. The power-law y index falls, or 
softens, with the temperature, see Fig. [161 An interesting feature 
appears to occur at the transition of the two power-law describ- 
ing the observed decrease of the temperature. The large time lag 
observed in the first episode and reported in section 6.1 has a 
clear explanation in the power-law behavior of the temperature 
and corresponding evolution of the photon index y. Figs. 1161171 



7.2. The radius of the emitting region 

We turn now to estimate an additional crucial parameter for the 
identification of the nature of the blackbody component: the ra- 
dius of the emitter rem- We have proved that the first episode is 
not an independent GRB, not a part of a GRB. We can there- 
fore provide the estimate of the radius of the emitter from non- 
relativistic considerations, just corrected for the cosmological 
redshift z. We have, in fact, that the temperature of the emitter 
Tern = Tobsi^ + z), and that the luminosity of the emitter, due to 
the blackbody emission, is 



L = 47rr,2^crrl = AtitI^ctTUI + z)\ 



(6) 



where Vem is the radius of the emitter and cr is the Stefan 
constant. From the luminosity distance definition, we also have 
that the observed flux ^obs is given by: 



l>obs 



47rZ)2 



Z)2 



We then obtain 

vl/2 



^ em 



9^obs 



D 



(8) 



The above ra dius differs from the radius Vph given in Eq. (1) 
of lRvde & Pe'erl (^009) and clearly obtained by interpreting the 
early evolution of GRB 970828 as belonging to the photospheric 
emission of a GRB and assuming a relativistic expansion with a 
Lorentz gamma factor F: 



rph = 



(1.06)(l + z)2 



(9) 



where *k = (^(pobsKcrT^^)^ ^ and the prefactor 1.06 arises from 
the dependence of rph on the angle of sight ('Pe'er,"2008). Typical 
values of rph are at least two orders of magnitude larger than our 
radius rem- We shall return on the analysis of GRB 970828 in a 
forthcoming paper. 

Assuming a standard cosmological model (Hq = 70 
km/s/Mpc, flm = 0.27 and Qa = 0.73) for the estimate of the 
luminosity distance D, and using the values for the observed flux 
(pobs and the temperature kTgbs, we have given in Fig. [181 the evo- 
lution of the radius of the surface emitting the blackbody rem as 
a function of time. 

Assuming an exponential evolution with time of the ra- 
dius in the comoving frame, we obtain from a fitting procedure 
the value S = 0.59 ± 0.1 1, well compatible with 6 = 0.5. We also 
notice a steeper behavior for the variation of the radius with time 
corresponding to the first 10 s, which corresponds to the emis- 
sion before the break of the double power-law behavior of the 
temperature. We estimate an average velocity of v = 4067 ±918 
km/s, R^ = 0.91, in these first 10 s of emission. In episode 1 the 
observations lead to a core of an initial radius of ~ 12000 km 
expanding in the early phase with a sharper initial velocity of ~ 
4000 km/s. The eff'ective Lorentz F factor is very low, F - 1 ~ 
10-^ 

8. Conclusions 

GRB 090618 is one of the closest (z = 0.54) and most ener- 
getic (Eiso = 2.9 X 10^^ ergs ) GRBs up to date. It has been 
observed simultaneously by the largest number of X and y ray 
telescopes: Fermi, Swift, AGILE, Konus-WIND, Suzaku-WAM 
and the CORONAS -PH0T0N-RT2. These circumstances have 
produced an unprecedented set of high quality data as well as the 
coverage of the instantaneous spectral properties and of the time 
variabilit y in luminosity of selected bandwidth of the source, see 
e.g. Figs. lll2[ In addition there is also the possibility of identify- 
ing an underlying supernova event from the optical observations 



(7) 
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Fig. 17. Evolution of the kT observed temperature of the black- 
body component and the corresponding evolution of the photon 
index of the power-law. The blue line in the upper panel corre- 
sponds to the fit of the time evolution of the temperature with a 
broken power-law function. It is evident a break time around 
lis after the trigger time, as obtained from the fitting procedure. 



in the light curve of well-defined bumps, as well as from the 
correspectiv e change in colou r after around 10 days from the 
main event (ICano et al.L l201lb . Unfortunately a spectroscopic 
confirmation of the presence of such supernova is lacking. We 
have restricted our attention in this paper to the sole X and y ray 
emission of the GRB, without addressing the possible supernova 
component. 

By applying our analysis within the fireshell scenario, see 
section 4, we have supported th at GRB 090618 is ac tually com- 
posed of two diff'erent episodes (iRuffini et al.Ll2010ah : episode 1, 
lasting from to 50 s and episode 2 from 50 s to 151 s after the 
trigger ti me. We have also illu strated the recent conclusions pre- 
sented in "Rufl&ni et al J (1201 ll) . that episode 1 cannot be either a 
GRB nor a part of a GRB, see section 5. By a time-resolved spec- 
tral analysis we have fitted the instantaneous spectra by a black- 
body plus an extra power-law component. The temperature of 
the blackbody appears to have a regular dependence with time, 
described by two power-law functions: a first power-law with 
decay index dikr = -0.33 ± 0.07 and the second one with hkr = 
-0.57 ± 0.11, see Section 7 . All these features follow precisely 
some of the results obtained by Felix Ryde and his collabora- 
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Fig. 18. Evolution of the radius of the first episode emitter, as 
given by Eq. dg]). 



tors (^ Rvde & Pe'eil|2009|) . where the authors analyzed selected 
temporal episodes in some GRBs observed by BATSE. 

We have also examined with particular attention, see section 
6, the radius rem of the blackbody emitter observed in the first 
episode, given by Eq. ([8]). We interpret the nature of this episode 
1 as originating froni what we have defined a proto-black hole, 
(Ru ffini et al.L 1201 Oah : the collapsing bare core leading to the 
black hole formation. Within this interpretation, the radius rem 
depends only on the observed energy flux of the blackbody com- 
ponent (pobs, the temperature kT as well as on the luminosity 
distance of the source D. We obtained a radius of the emitting 
region smoothly varying between ~ 12000 and 70000 km, see 
Fig-im Other interpretations associating the origin of this early 
emission to the GRB main event (Pe'er et al., 2007) lead to a dif- 
ferent definition for the radius of the blackbody emitter, which 
results to be larger than our radius by at least two orders of mag- 
nitude. We are planning a systematic search for other systems 
presenting these particular features. 

Episode 2 is identified as a canonical long GRB which orig- 
inates from the black hole formation process and lasts in arrival 
time from 50s to 151 s after the trigger time. The good qual- 
ity of data allowed us to search for the P-GRB signature in the 
early emission of the episode 2. From a detailed analysis we find 
that the first 4 s of episode 2 are in good agreement with the 
theoretically predicted P-GRB emission, see section 5.2. The 
observed spectrum integrated over these 4 s is well fitted by a 
blackbody with an extra power-law component, where this latter 
component is mainly due to the early emission of the extended- 
afterglow, see Fig. [H From the temperature observed in the P- 
GRB, kTpGRB = 29.22 ± 2.21, and the E^^f energy of the sec- 
ond episode, which we assumed equal to the isotropic equiva- 
lent energy of this episode, E^^f = 2.49 x 10^^ ergs, we ob- 
tained the value of the baryon load of the GRB, see also Fig. 
H 5 = (1.98 ± 0.15) X 10"^ and a consequent Lorentz F fac- 
tor at the transparency of Fo = 495 ± 40. We have been able 
to simulate the temporal and the spectral emission of the sec- 
ond episode, as seen by the Fermi-GBM instrument (8 keV - 10 
MeV). As we have shown in Fig. [121 our simulation succeeds 
in fitting the light curves as well as the spectral energy distribu- 
tion emitted in the first main spike of the second episode. The 
residual emission of the last spikes is reasonably fitted, taking 
into due account the difl&culties in integrating the equations of 
motion, which after the first interactions of the fireshell with the 
CBM become hardly predictable. The energetic of the simula- 
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tion is fulfilled and we find that the emission is due to blobs of 
matter in the CBM with typical dimensions of r/^/ = 10^^ cm 
and average density contrast dnln^l particles/cm^ in an over- 
all average density of 1 particle/cm^. We need to find additional 
cases of such phenomena to augment our statistic and improve 
its comprehension. 

Particularly relevant are the first two-dimensional hydrody- 
namical simulations of the progenitor evolution of a 23Mq star 
close to core-collapse, leading to a naked core, as sh own in the 
recent work of Arnett and Meakin (lArnett & Meakin . 2011). In 
that work, pronounced asymmetries and strong dynamical inter- 
actions between burning shells are seen: the dynamical behav- 
ior proceeds to large amplitudes, enlarging deviations from the 
spherical symmetry in the burning shells. It is of clear interest to 
find a possible connection between the proto black hole concept, 
introduced in this work, with the Arnett and Meakin results: to 
compare the radius, the temperature and the dynamics of the core 
we have found in the present work with the naked core obtained 
by Arnett and Meakin from the thermonuclear evolution of the 
progenitor star. Particularly relevant is the presence, during this 
phase of collapse, of strong waves, originated in the mixing of 
the diff'erent element' shells. Such waves should become com- 
pressional, as they propagate inward, but they should also dis- 
sipate in non-convective regions, causing heating and slow mix- 
ing in these regions of the star. Since the wave heating is faster 
than radiative diff'usion (which is very slow), an expansion phase 
of the boundary layers will occur, while the iron (Fe) core will 
contract (^Arnett & Meakinl |2Q11|) . There is also the interesting 
possibility that the CBM clouds observed in GRBs be related to 
the vigorous dynamics in violent activity of matter ejected in the 
evolution of the original massive star, well before the formation 
of the naked core (Arnett D., private communication). 

It is appropriate to emphasize that these results have no re- 
lation with the study of precursors in GRBs done in the current 
literature (see e.g. Burlo n et al.L 12008L and references therein). 
Episode 1 and episode 2 are not temporally separated by a qui- 
escent time. The spectral feature of episode 1 and episode 2 are 
strikingly diff'erent and, moreover, the episode 1 is very ener- 
getic, which is quite unusual for a typical precursor event. We 
finally conclude that for the first time we witness the process of 
formation of the black hole from the phases just preceding the 
gravitational collapse all the way up to the GRB emission. 

There is now evidence that the Proto Black Hole forma- 
tion has been observed also in other GRB sources. After the 
submission of this article a second example has been found in 
GRB 101023, then and a paper about this source was submit- 
ted on November 4th 2011 and then published on February 1st 
2012 (Penacchioni et al.","2012). There, extremely novel consid- 
erations in the structure of the late phase of the emission in X- 
ray at times larger than 200 s have been presented in favour of a 
standard signature in thes e sources (see also the considerations 
made in lPage et"ail |201 ik . The possible use of this new family 
of GRBs as distance indicators is being considered. 
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